Dynamo in Protostar 

MAHENDRA K. VERMA^' *, BIDYA BINAY KARAK''' t, ROHIT 
KUMAR^í 

Department of Physics, Indián Institute of Teciinology, Kanpur 208016, India 
''Department of Physics, Indián Institute of Science, Bangalore 560012, India 

Abstract. In tliis páper, we estimate tfie magnetic Reynolds number of a typical protostar before 
C^ ■ and after deuterium buming, and claim for the existence of dynamo process in both the phases, 

because the magnetic Reynolds number of the protostar far exceeds the critical magnetic Reynolds 
^^ , number for dynamo action. Using the equipartition of kinetic and magnetic energies, we estimate 

the steady-state magnetic field of the protostar to be of the order of kilo-gauss, which is in good 
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5-H ' agreement with observations. 
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1. Introduction 
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rS^ ' Magnetic field is omnipresent in the universe. It is found in many stars, planets, galaxies, 

Mh. interstellar medium, etc, and it is believed to induce inflow of matter or accretion in 

Q I astrophysical disc and immersed cloud, coUimated jets and outflows, magnetic breaking, 

and protostellar winds. The magnetic field plays a crucial role in all stages of the evolution 
C/3 ' of a protostar. In this páper we will address an important question related to protostars: 

when does the magnetic field appears in protostars, and what is the magnitude of the 
magnetic field, if present? 

Presence of X-ray from the protostar regions suggests a strong magnetic field in its 
centrál region Ju |2|] . Also, the detection of polarized synchrotron emission arising from 
^^ . protostar jets [3] and the measurements from Zeeman broadening of photospheric lineš 

(yQ ' 14-7J provide a strong support for the existence of the magnetic field in the protostars. 

^O ■ Several questions on protostars are: whether the magnetic field arises from the molecular 

T- 1- I cloud of the protostar formation region, or due to šelf induction or dynamo mechanism? 

f^ . How is the magnetic field sustains itself and compensates its removal due to the outflows, 

CO ' or its dissipation as ohmic heating? Does the dynamo mechanism, if any, stops after the 

formation of the protostars? 

Machida et al. [8J and Sur et al. (|9|] argue that the dynamo mechanism can amplify the 

i!^ . initial magnetic field of the stár forming regions. On the other hand, Tan and Blackman 

'V^ I OlOll suggest that dynamo amplification of the primordial magnetic field is possible in 

C^ . protostellar disc, and they propose that the induced magnetic field is helical. In these 

work, the generation mechanism for the magnetic field involves several mechanisms, e.g., 

convection, rotation, etc. However, there is no substantive statement in the literatuře on 

the necessity of dynamo process in protostars. In the present páper we claim that the 

protostars are magnetic because its magnetic Reynolds number to tiigger the dynamo is 

far above the critical magnetic Reynolds number observed in laboratory expeiiments and 

numerical simulations. The equipartition of kinetic and magnetic energies yields magnetic 
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field of the protostar to be of the order of kilo-gauss, which is in good agreement to those 
observed in nature. We present these arguments in the next two sections of the páper, and 
conclude in the last section. 

2. Critical values of magnetic Reynolds number for dynamo 

The fluid velocity u and the magnetic field B in a dynamo mechanism are govemed by 
the magnetohydrodynamic (MHD) equations: 

p{dtu + (u • V)u} = -\7p + (J X B) + //V^u + F, (1) 

dtB = Vx{uxB)+riW^B, (2) 

V • B = O, (3) 

where u, J, F, p, v, and 77 represent the velocity field, current density, external forcing, 
hydrodynamic pressure, kinematic viscosity, and magnetic diffusivity, respectively. Notě 
that J = c(V X B)/47r. The density of the fluid is governedby the continuity equation. 



though some computations assume it to be a constant 11 1111 . Notě that only the velocity 
field is being forced in the above equation. 

Dynamo mechanism is said to occur when the magnetic energy reaches a finite value 
asymptotically, that is as í — ?► 00. The magnetic field gets energy from the velocity field 
through the V x (u x B) term, and loses energy via Joule dissipation. For magnetic 
energy to grow and reach a steady statě, a comparison of these two terms suggests that 

(V x (u X B)) > T^V^B (4) 

or, the magnetic Reynolds number 

Rm = — > 1 (5) 

The minimum magnetic Reynolds number, called the critical magnetic Reynolds number 
Rmc, could be larger than one due to geometrical consideration or other factors like mag- 
netic Prandtl number (Pm, defined as v/rj), rotation, etc. Yet, we can argue that Rmc is 
of the order of unity for all dynamo processes; this limit has been consistently observed 
in a number of laboratory experiments and numerical simulations. 



Among laboratory experiments, Gailities (Riga experiment) et al. Ill2ll . Stie ghtz and 



Miiller (Karlsruhe experiment) 111311 . and Monchaux et al. (VKS experiment) 111411 ob- 
served self-generated magnetic fields in their experimental setups. Liquid sodium, whose 
magnetic Prandtl number is approximately 10^^, was ušed as the operating fluid in the 
aforementioned experiments. The critical magnetic Reynolds number for the above set of 
experiments was greater than 10; specifically, the Rmc ~ 30 for that Monchaux et al.'s 
experiment, commonly referred to as VKS experiment 111411 . 

A large number of high-resolution simulations háve been performed to study the dy- 
namo transition. Most researchers ušed pseudo-spectral method to solve the MHD equa- 
tions for various forcing in a box geometry. The range of magnetic Prandtl number ušed 
so far is from 10^^ to 10^. Notě that numerical simulations of very large Pm or very 
small Pm dynamo are difficult due to their requirements of high resolutions. Notably, 
Schekochihin et al. Í15[ll6ll . Ponty et al. Í17I1 . and Iskakov et al. \ 18] simulated dynamo 



and studied variation of Rmc as a function of magnetic Prandtl number They observed 
dynamo for both small and large Prandtl numbers, with the range of Rmc between 10 
and 500. The Rmc for lower Prandtl number tends to higher in most of the simulations. 
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Recently, Yadav et al. Í19ll20tl studied Taylor-Green dynamo and observed the Rnic to be 
around 10. In Fig. 1 we plot, Rnic for the aforementioned experiments and numerical sim- 
ulations. We can conclude from these results and the aforementioned phenomenological 
arguments relating the nonlinear term to the dissipative term that Rmc varies with geome- 
try, forcing, and Prandtl number, but it remains bounded between 1 and 500. Hence, if the 
protostars háve Rm > Rmc, then dynamo must be active in the protostar. This is what 
we would estimate in the next section. 
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Figuře 1. Plot of critical magnetic Reynolds number (Riiic) vs. magnetic Prandtl 
number (Pm). The plot exhibits data of VKS experiment I14ll . Riga experiment U2ll. 
Karisruhe experiment Il3ll , and numerical data of Ponty et al. DNS |17], Yadav et 
al. DNS [il,|23], and Schekochihin etal. DNS L15.,i60. 



3. Magnetic Reynolds number and dynamo action in the protostar 



In this section we will discuss the dynamo process in protostar before and after deuterium 
burning. First, we focus on a typical protostar before deuterium burning. Temperature of 
the protostar at this stage may be around 10^ K, its rádius around 100 solar-radii, and its 
mass around two solar mass. The kinematic viscosity and the resistivity of the gas can 
be estimated using Spitzer formula (see e.g., Schekochihin et al. Í2lll or Choudhuri MM) 
that provides 



v - 2.21 X 10-i^r^/V(4/9), 

7? ~ An^/^m}J^e^c^l{2 x {2kBTf'^ x 0.6), 



(6) 
(7) 



where p is the density of the medium, T is the temperature in Kelvin, fc^ is the Boltz- 
mann's constant, TTig and e are the mass and charge of an electron, and c is the speed of 
light. Notě that the above formulas are in CGS units. We substitute the aforementioned 
parameters, T = 10^ K and the average density p sa 3 x 10^^ gm/cm^^, in the above 
formulas. Also, we také the large-length scale L of the protostar to be 10^° cm, and the 
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large-scale velocity U to be one-tenth of the sound speed. Consequently, 

1^ K 2 cmVs, (8) 

íí w 3 X 10^ cmVs, (9) 

U « 10^ cm/s, (10) 

Rm ft! 5 X 10^, (11) 

Pm « 10"^ (12) 

The magnetic Reynolds number Rm is much greater than the estimated critical mag- 
netic Reynolds number (10 — 500) using the laboratory experimental and numerical data. 
Hence we expect dynamo to be active in typical protostars before the deuterium buming 
stage. The background galactic magnetic field or the magnetic field amplified during the 
formation of protostars could act as the seed magnetic field for this process. 

The saturated or the steady-state magnetic field of the proto stár can be estimated easily 
using the equipartition of kinetic and magnetic energy, i. e., B"^ /Att ~ pU^, which yields 
the average magnetic field at this stage to be around 700 Gauss. Notě that approximate 
equipartition of kinetic and magnetic energy in the steady statě has been observed in 
large number of systems (e.g., solar wind) and numerical simulations, hence it is a robust 
assumption for estimating the steady-state magnetic field. Numerical simulations also 
reveal that saturation of the magnetic field takés around 10 eddy turnover time. Hence we 
expect the time taken for the magnetic field of the protostar to reach a steady statě would 
be approximately lúL/U , which is approximately 12 terrestrial days. 

Now let US study the possibility of dynamo process in the protostar after deuterium 
burning. At this stage, the temperature of a typical protostar is approximately 10^ K, and 
its rádius around two solar-radius. Here, we také the large-length scale L to be one tenth 
of the solar rádius, and time period of the large eddies to be the rotation time-period of 
the protostar, which is approximately 20 terrestrial days. Hence, 

(13) 
(14) 
(15) 
(16) 
(17) 

Again, since Rm > Rmc, we expect dynamo to be active in the protostar after the deu- 
terium burning phase. The magnetic field generated in the pre deuterium-burning stage 
would act as the seed magnetic field. Equipartition of the kinetic and magnetic energy 
yields the average magnetic field to be of the order of 2 x 10^ Gauss, and it would také 
around 200 terrestrial days for the magnetic field to reach saturation. Notě that our simple 
estimate of the protostar magnetic field is in generál agreement with many observational 
measurements of Zeeman broadening |4-7], which reveal the existence of magnetic field 
of kilogauss strength. The high field thus generated could produce X-ray activity, jets, 
outflows, and other magnetic effects. We remark that the above estimates are in good 
agreement with the observed parameters of the Sun. 

4. Conclusions 

In this páper, we estimate the magnetic Reynolds number Rm of a typical protostar before 
and after deuterium burring. We show that the Rm in both the phases is far greater than the 
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critical magnetic Reynolds number Rnic, which could be between 10 and 500, an estimate 
based on dimensional analysis, laboratory experiments, and numerical simulations. Thus, 
we claim that the dynamo mechanism is active in protostar in both the phases. We also 
estimate the steady-state magnetic field of the protostar in these two phases using the 
equipartition of kinetic and magnetic energy, and find them to be of the order of a kilo- 
gauss, which are in very good agreement with several astronomical observations |4-7]. 

Our arguments on protostar dynamo is quite robust, and it is independent of mechanism 
invoked for the process, e.g., convection, rotation, a — cj, etc. It must be however kept in 
mind that the actual magnetic field of the stár would depend on the details of the dynamo 
process. Our arguments provide an argument in favour of the dynamo mechanism, as well 
as an estimate of the strength of the magnetic field, in protostar. 
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